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1.  Introduction  and  problem  formulation 

Using  the  Laplace  transform,  we  (Jordan  and  Puri,  1999)  derived 
the  exact  solution  to  Stokes’  first  problem  for  a  dipolar  fluid  (see 
Cowin  (1974),  Straughan  (1987)  and  the  references  therein)  that 
is  valid  for  arbitrary  values  of  the  material  parameters  d(>0)  and 
/(>0).  The  present  Note  is  devoted  to  extending  this  earlier  work 
of  ours  to  the  case  of  micropolar  fluids,  the  theory  of  which  was 
formulated  by  Eringen  (1966).  In  doing  so,  we  also  extend  the  re¬ 
sults  of  Jordan  et  al.  (2000)  on  Laplace  inverses  involving  doubly- 
nested  square  roots  with  one  branch  point  to  a  more  general  class 
of  such  inverses  possessing  two  branch  points. 

To  this  end,  we  observe  that  by  replacing  the  straight  channel  of 
Kirwan  and  Newman  (1972), 2  who  considered  the  unsteady  plane 
Couette  flow  of  an  incompressible  micropolar  fluid,  with  the  half¬ 
space  y  >  0,  and  setting  the  pressure  gradient  term  in  their  equation 
of  motion  to  zero,  one  possible  way  to  formulation  Stokes’  first  prob¬ 
lem  for  a  micropolar  fluid  is  as  the  IBVP 

pJun  +  2kut  -  k(2v  +  K)uyy  -  [}>  +  J(p  +  k)]uwt 

+  y(V  +  K)Uyyyy  =  0,  (j/,  t)  6  (0,  00)  X  (0,  00) ,  (1.1) 

u(0,t)  =  U00(t),  u(oo,t)  =  0,  uyy(0,t)  =  Ms,  uyy(oo,t)  =  0  (t  >  0) , 

(1.2) 
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2  The  mass  density,  p,  appears  to  have  been  inadvertently  omitted  from  Eqs.  (7)-(9) 

of  Kirwan  and  Newman  (1972);  see  Devakar  and  Iyengar  (2009),  wherein  the  first  two 
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u(y,0)  =  ut(y,0)  =  0  (y  >  0).  (1.3) 

In  IBVP  (1),  the  velocity  vector  has  the  form  v=  (u(y,  f),  0,  0), 
where  the  incompressibility  condition  V  ■  v  =  0  is  identically  satis¬ 
fied;  0(  )  denotes  the  Heaviside  unit  step  function;  y,  k,  and  the 
dynamic  viscosity  ft  are  positive  constants;  here,  v  =  pip  denotes 
the  kinematic  viscosity;  U0  is  a  nonzero  constant;  Ms  is  taken  as 
constant;  and  we  have  set  k:  =  k/p  for  convenience. 

Remark  1.  The  boundary  conditions  (BC)s  assumed  in  Eq.  (1.2), 
which  in  the  present  context  are  equivalent  to  BCs  of  “type  A”  in 
dipolar  fluid  theory  (Straughan,  1987),  were  suggested  by  Kirwan 
and  Newman  (1972).  For  more  on  the  question  of  what  BCs  are 
appropriate  in  the  micropolar  context,  see  also  Eringen  (1966)  and 
Cowin  (1974). 


2.  Exact  solution  via  the  Laplace  transform 


To  reduce  the  number  of  coefficients,  we  introduce  the  follow¬ 
ing  nondimensional  quantities:  u  =  u/U0,  y  =  y/L0,  t  =  t/T0,  where 


Lo  = 


y  +  Jr(2p  +  k) 


and 


21 

t  ZLo 

i  o  =-. 


2k  ‘ u  2v  +  K  ’ 

and  recast  IBVP  (1)  in  the  simpler  form 

hutt  +  Ut-  liyy  -  Uyyt  +  PUyyyy  =  0,  (}/,  t)  6  (0,  (X))  X  (0,  OO), 

U(O,t)  =  0(t),  U(OO,t)=0,  Uyy(O.t)  =M,  Uyy(  OO  ,  t)  =  0 

u(y,0)  =  ut(y,0)  =  0  (y  >  0), 


(2) 


(3.1) 

(t>  0), 

(3.2) 

(3.3) 
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where  M  is  the  nondimensional  form  of  Ms  and  all  primes  have  been 
suppressed  but  remain  understood.  Furthermore,  t  >  0  and  h  g 
(0,1/2)  are  given  by 


1  y(fi  +  k)  S(2n  +  k) 

Lo  V  k(2j.i+  k)’  4 kLl  ' 


(4) 


On  applying  the  temporal  Laplace  transform,  se[-\,  to  Eqs.  (3.1) 
and  (3.2),  using  the  initial  conditions  given  in  Eq.  (3.3),  and  then 
solving  the  resulting  subsidiary  equation,  the  transform  domain 
solution  is  found  to  be 


u(y,s)  = 


t2  {e~r2yr\  -  e~r'yr2  -  M{e ^  -  e^) 


sJs2(  1  -  4hr)  +  2s(l  -  2t)  +  1 


(5) 


where  s  is  the  transform  parameter,  a  bar  over  a  quantity  denotes 
the  image  of  that  quantity  in  the  Laplace  transform  domain  (e.g., 

u(y,s)  :=  SC[u(y,t)]),  and 


ri,2=j 


s  +  1  ±  \J(s+  l)2  -  4 (2s(hs  +  1) 


(6) 


To  simplify  the  task  of  inversion,  it  is  convenient  to  recast  Eq. 
(5)  as 


u(y,s)  =  Uj(y,s)  -  (2 £2M  -  1  )u2(y,s)  +  u3(y,s), 

where 

e-w  _|_  e-ny  _  e-r2y  _  e-n y 


(7) 


u  1 


2s 


u  2  - 


2s>J(s+\)2  -  At2 s(hs  +  1)’ 


u3  =su2. 


(8) 


Since  the  inverses  we  require  do  not  appear  to  be  given  in  any  of 
the  standard  tables,  we  turn  to  the  complex  inversion  formula  (see, 
e.g.,  Duffy,  2004)  and  observe  that  the  inverses  of  the  terms  in  Eq. 
(7)  are  given  by 


Ui,2,3(y,t)  lim  /  estuli2,3(y,s)ds  (t  >  0), 

zm  t-kjo  ja_iT 


(9) 


where  Uj ,2,3(y,t)  =  0  for  t  <  0  and  the  arbitrary  real  number  cr  >  0  lies 
to  the  right  of  all  singularities.  The  first  step  in  evaluating  these 
integrals  is  determining  the  singularities  of  their  integrands.  Omit¬ 
ting  the  details,  it  can  be  shown  that  the  only  non-removable  singu¬ 
larities  of  u  12  consist  of  a  simple  pole  at  s  =  0  and  branch  points  at 
s  =  { — l//i,  0}  while  those  of  u3  consist  of  only  the  two  branch  points. 
Thus,  by  Cauchy’s  theorem  we  have 


jfestu,,2,3(y,s)ds  =  0. 


(10) 


where  the  Bromwich  contour  C  (see  Fig.  1)  is  traversed  in  the  coun¬ 
terclockwise  direction. 

Following  standard  techniques  of  contour  integration  (see,  e.g., 
Duffy,  2004),  it  can  be  established  that,  for  all  h  e  (0,1/2)  and  t  >  0, 

3  exp[-t(j/+  l//i)]sin[y<»+(f/)]cos[y(»_(i/)]dr/ 


n 


ri  +  1/h 

exp[- t(r\  + 1  /h)]  sintyco+(jj)] cos(yeo_ (i/)]df/ 


i/  +  l /h 

exp[-t(i;  +  l//i)]sin[y(n+(r/)]cosh()/co_(!7)]d)/ 
rj+\/h 


ey(4hy'/2, 
t>(4  Hyy2, 


MW)  - m  Id  ■ - e->«)  •  r  e-'smLv^cHd,, 

[  7°  i]\J (rj-\)2 +4r]l2(\ -hr\)^ 


m 


’  r°°  exp [— 1(»7  + 1  /h)]  sin[yoj_()j)]  cos[yco+(r/)]d)/ 
Jo 

L 


(tl  +  \/h)^[{i]+\/h)-\\2-4ril2(hri  +  \) 
exp[-t(>/+ 1  /h)]sin[y(B_(r;)]cos[yft)+(r;)]df; 

('7  +  1  /h) /h)  - 1}2  -4ijt2  (htj  +  1) 
exp[-t(/;+ 1  /h)]sinh[ytu_(t/)]  cos[yco+  (»7)]  d/y 
’r  (r\+\/h)\j\[{r]+\/h)-\f  -4rji2(ht]+\)\ 


d 


ty(4h) 


e>(4h) 


-1/2 


-1/2 


(12) 


rVh 


(  Jo  \m-l)  +4;?f2(l-fH7)J 


m 


'  r  exp[-t(»;  + 1  /7i)]  sin[yoj_(/7)]  cos[yco+(j/)]dj/ 
Jo 

L 


\J[(t]+\/h)  - 1]2  -  4r]l2(ht]  + 1 ) 
i'  exp[-f(»7+ 1  /h)]sin\ya>_(ri)]cos\yw+(ri)]di] 
\Z[('/  +  1/1')  -  l]2  -4j/f2(/u/  + 1) 
exp[-t(r/+ 1  /h)]sinh[yco_(^)]  cos\yco+(rj)]di] 


t<{4h) 


-1/2 


d 


\J\[(n  +  l/h)- 1}2 - 4i]l2 (hi)  +  ))\ 


where 


i]  -  1  +  \J (rj  -  l)2  +  4h£2r/(l/h  -  i]) 


e>(4hy'/2-, 

(13) 

(14) 

(15) 

(16) 


C0±^  =  Tt\j  (f/  +  V^)  -  1  2£\Jht](rj  +  J/h) 

,  _  1  -  h(l  +  2f2  -  2lVl2yh-\) 

11  ~  h(4h(2  -  1)  ’ 

and  we  observe  that  rf  >  0  if  and  only  if  t  >  1  /V4h. 

3.  Special/limiting  cases  and  other  applications 


Remark  2.  The  steady-state  solution,  u^y),  is  defined  as 
Uoo(y):  =  limt_>00u(y,t)  and  is  given  by  u^y):  =  1  -  M£2[l  - 
exp ( — y/^)],  which  is  identical  in  form  to  that  of  the  dipolar  fluid 
case  (see  Jordan  and  Puri,  1999,  Eq.  (4.3)). 

Remark  3.  If  we  take  h  =  1  - 12,  implying  t2  e  Q,l),  then  u1-2,3 
simplify  to  such  an  extent  that  they  can  be  easily  inverted  using, 
e.g.,  the  tables  of  Laplace  inverses  given  in  (Erdelyi,  1954,  Chap. 
V).  Thus,  for  this  special  case  u(y,  t)  =  <W(y,  t),  where  *  is  defined  as 


at 


My  r;  :=fl(()|erfc^+_)  2  e“ 

_|_e-y\/5erfc^_Z__  _  ecyVa+5erfc^^^+  (a  +  b)tj 

_e-cyvraerfc^_v/(a  +  b)^ 

-  Mt2  erf cf ^  eat  ey'/5erfcf -^-= +  Vat] 

[  \2  Vi)  2  [  V2  Vi  ) 

+  e  yy/iedc^^j=  -  VatJ  -  eq''/S+5erfc^^=  +  ^/(a  +  b)tj 
_  e-o/i/a+Eerfc  -  \J(a  +  fyt'J  -  ^  e^erfc  ^ ^ =  +  Vbt'j 

i  ^  cy'+erfc  ^  ^  ,/m^)  l  ( i  vi 
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Im(5) 


Re(s) 


Fig.  1.  Bromwich  contour,  C,  used  to  evaluate  the  integrals  in  Eq.  (10).  Here,  e  and  e  are  small  positive  quantities. 


Here,  erfc  (•)  denotes  the  complementary  error  function  and  the 
positive  constants  a,  b ,  and  c  are  defined  as  a  =  (2£2  -  l)-1,  b  = 
(1  -  f2)-1,  and  c  = 

Remark  4.  Setting  £  =  l2  and  letting  h  ->  0  reduces  u  to  the  l\  =  \ 
special  case  of  (Jordan  and  Puri,  1999,  Eq.  (3.4)),  i.e.,  the 
corresponding  solution  for  the  case  of  a  dipolar3  fluid.  On  the 
other  hand,  letting  £  -»  0  reduces  u  to  the  solution  of  Stokes’  first 
problem  for  an  Oldroyd-B  liquid;  see  pp.  354-355  of  Duffy  (2004), 
and  also  Christov  and  Jordan  (2009). 

Remark  5.  It  is  noteworthy  that  Eq.  (1.1),  which  is  identical  in 
form  to  the  equation  of  motion  for  this  problem  involving  a  micro- 
structure-laden  fluid  of  the  type  described  by  the  theory  of  Kline 
and  Allen  (1970),  also  arises  in  the  study  of  incompressible  flow 
in  fractured  or  fissured  media  (Guenther  and  Lee,  1996).  Eq.  (1.1) 
can  also  be  regarded  as  the  governing  equation  for  a  particular  type 
of  two-species  reaction-diffusion  system. 
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3  In  the  context  of  the  present  problem,  dipolar  fluid  theory  and  Green-Naghdi 
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see,  e.g.,  Quintanilla  and  Straughan  (2005)  and  the  references  therein. 
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